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ABSTRACT 

High cadence, multi-wavelength observations and simulations are employed for the analy- 
sis of solar photospheric magnetic bright points (MBPs) in the quiet Sun. The observations 
were obtained with the Rapid Oscillations in the Solar Atmosphere {ROSA) imager and the 
Interferometric BIdimensional Spectrometer {IBIS) at the Dunn Solar Telescope. Our anal- 
ysis reveals that photospheric MBPs have an average transverse velocity of approximately 
1 kms~^, whereas their chromospheric counterparts have a slightly higher average velocity 
of 1.4 kms~^. Additionally, chromospheric MBPs were found to be around 63% larger than 
the equivalent photospheric MBPs. These velocity values were compared with the output of 
numerical simulations generated using the MURaM code. The simulated results were similar, 
but slightly elevated, when compared to the observed data. An average velocity of 1 .3 km s^^ 
was found in the simulated G-band images and an average of 1.8 kms~^ seen in the velocity 
domain at a height of 500 km above the continuum formation layer Delays in the change of 
velocities were also analysed. Average delays of s between layers of the simulated data set 
were established and values of ^29 s observed between G-band and Ca II K ROSA observa- 
tions. The delays in the simulations are likely to be the result of oblique granular shock waves, 
whereas those found in the observations are possibly the result of a semi-rigid flux tube. 
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1 INTRODUCTION 

Magnetic bright points (MBPs ) are ubiquitous in the solar pho- 
tosphere l lDunn & ZirkeJ[r973l) . They are thought to be the foot- 
points of magnetic flux tubes in the in ter-granular lanes with 
field strengths of the order of a kilogauss dStenflollTgsl :! Solanki I 
Il993l) . Their brightness is due to the reduced pressure within the 
flux tube allowing the observer to view a deeper, hotter region 
of the photosphere, and is also due to the heating of the plasma 
within the flux tube by material surrounding its walls. The higher 
temperature reduces the abundance of the CH mol ecule, making 
MBPs appear brighter in G-band (4305 A) images l lSteiner et al. I 
Shelyag et al. 2004). G-band observations therefore provide 
a good base for analysing the photospheric properties of these con- 
centrated magnetic structures. 

The importance of MBP studies was reviewed in 
Ide Wiin et al.T (12009'). 
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transmitted to the corona and lead to he ating. This has been a 
recent topic of several 3-D MHD models 1 van Ballegooiien et al. I 



l201ll : lAsgari -Targhi & van Ballegooiien 1l2012h investigating these 
effects. 

Recent MBP studies have focussed on their size, velocit 



distribution s 1 Sanchez Almeida et al. Il2004|: Crockett et al 



Kevs et al] I2OI1I) and intensity peaks jLanghans et al. I 20021) . 
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Bonet et al. I ( l2008 h have identified convective downdrafts and vor- 



Briefly, studies of 
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MBP 



motions are 



important as these may drive Alfven waves in magnetic field lines 
within flux tubes, where the energy carried by the waves can be 
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tex motions at the location of MBPs, phenomena that have also 
been reproduced in simulations of radiative MHD jShelvag et al. I 
l2011al : IShelvag et al.ll2011bll2012h . 

A study of MBPs in the NalDi line jjess et al.ir2010al) has 
shown evidence for strong downdrafts with velocities as high at 
7 km s~^ at MBP centres. It was also noted that field strengths af- 
fected the expansion of the flux tubes, with weak fields leading to 
an expansion of ~76% whereas higher field strengths led to an ex- 
pansion of ~44%. 

iKevs et alj ( l201lh studied the transverse velocities of over 
6000 MBPs in both observed and simulated G-band images. Their 
study revealed an average velocity of 1 kms^^ and a maximum 
value of 7 kms^^. Radiative MHD simulations with an average 
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Figure 1. Sample ROSA and IBIS data employed in this study, (a) G-band image (b) Co-spatial and co-temporal image in Call K. (c) and (d) Co-spatial IBIS 
images in the wing and the core of the Na I Di line, respectively. MBPs are abundant in (a). However, higher in the atmosphere their boundaries become 
more ambiguous due to reverse granulation. (e) An expanded view of the 4 data sets in the region outlined by the white box in (a), (i), (ii), (iii) and (iv) show 
expanded regions in G-band, Ca II K and the wing and core of the Na I Di line, respectively. A small sample of MBPs are contoured red in the 4 data sets 
showing both groups of MBPs and isolated MBPs. The different numbers and sizes of MBPs is evident across the 4 data sets. 
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field strengths of 200 G and 400 G revealed similar velocity distri- 
buti ons, the latter showing the b est agreement with observations. 

lUlmschneider et alj ( Il99lh have shown that the horizontal 
velocity of solar magnetic flux tubes increases as a function of 
height. The amplitude of the horizontal velocity component scales 
as p'^^^, where p is the density, and the authors note that the trans- 
verse wave flux is only approximately conserved. In the regime 
of strict flux conservation, a horizontal velocity component of 
1.5 kms~^ was expected at a height of 340 km, yet a velocity of 
only 0.8 kms~^ w as observed. This is very similar to the values of 
iKevs et al.l lioTH) . 

MBPs have average diameters of around 135 km 
jSanchez Almeida et al. I |2004|) and lifetimes of approximately 
90 s jKevs et alj 2011 ). and hence tracking these throughout the 
solar atmosphere requires high spatial and temporal resolution 
imaging. The accurate co-alignment of the images is vital when 
multi-wavelength studies are considered. In this paper we use 
simultaneous multi-wavelength observations to study MBPs at 
various heights in the solar atmosphere and improved BP tracking 
algorithms over our previous studies (Crockett et al. 2010). We 
present multi-wavelength observations and numerical simulations 
of photospheric MBPs in co-temporal and co-spatial G-band, 
Ca II K and Na I Di images. The observations and numerical 
simulations used are described in Section |2] Analysis of these 
results as well as their implications are presented in Section[3] with 
concluding remarks given in Section|4l 



2 OBSERVATIONS AND NUMERICAL SIMULATIONS 

The imaging data employed in this study were obtai ned using 
the Ra pid Oscillations in the Solar Atmosphere ( ROSA. ' Sess et al.l 
l2010bh instrument, at the Dunn Solar Telescope (DST), in New 
Mexico, USA. Observations were obtained during a period of ex- 
cellent seeing on 2009 May 28, using a 9.2 A wide filter centred at 
4305 A (G-band) and a 1 A wide filter centred at 3933.7 A (Ca II K 
core). We observed a 69" x 69" quiet Sun region at disk centre for 
~50 minutes, using a common plate scale of 0" .069 pixen ^ for the 
ROSA camer as. The images were reconstructed using the Speckle 
algorithm of IWoger et alj ( |2008|) . while image de-stretching was 
performed using a 40 x 40 grid (equating to a «l".7 separation 
between spatial samples; Jess et al. 2008). G-band images were 
taken at a raw cadence of 0.033 s, while after reconstruction the ca- 
dence was reduced to 0.528 s. Reconstructed G-band images were 
then binned into consecutive groups of four to improve the signal- 
to-noise and reduce the overall volume of the dataset, providing a 
final image cadence of 2.1 s. The Ca II K images had a raw ca- 
dence of 0.26 s, and a post-speckle reconstructed cadence of 4.2 s. 
Subsequently, the two data sets were spatially co-aligned. 

The Interferometric BIdimensional Spectrometer (IBIS; 

ICava llini 2006) was employed in conjunction with ROSA, cover- 
ing the Na I Di absorption line at 5895.94 A. At a spatial sam- 
pling of 0".097 pixeP^, the near square field-of-view of ROSA 
was contained within the circular aperture of IBIS. Nine wavelength 
steps were utilised with IBIS, with 10 images per wavelength step 
to aid image reconstruction, resulting in a complete scan cadence 
of 39.7 s. A blues hift correction was required due to the classical 
etalon mountings j Cauzzi et 'Zl l2008f) . Images of the Na I Di core 
used for MBP studies are true intensity maps, created by establish- 
ing the line-profile minimum at each pixel. By displaying Doppler- 
compensated line-centre intensities, rather than rest-wavelength in- 
tensities, brightness variations throughout the image are more in- 



dicative of the source function than of the velocities present in the 
line-forming region ( Jess et al. 2010a). A sample of ROSA and IBIS 
images is shown in Figure [T] 

The radiative m agneto-hydrodynamic (MHD) code MuRAM 
dVogler et aLll2005h . has been used to produce simulated G-band 
images, and their corresponding velocity domain. This code solves 
large-eddy radiative three-dimensional MHD equations on a Carte- 
sian grid, and employs a fourth-order Runge-Kutta scheme to ad- 
vance the numerical solution in time. The numerical domain has a 
physical size of 12x 12 Mm^ in the horizontal direction, 1.4 Mm 
in the vertical direction, and is resolved by 480x480x100 grid 
cells, respectively. Our starting point for the simulations is a well- 
developed non-magnetic (B — 0) snapshot of photospheric con- 
vection taken approximately 2000 s (about 8 convective turnover 
timescales) from the initial plane-parallel model. A uniform mag- 
netic field of 200 G was introduced at this stage, and a sequence of 
829 snapshots recorded, each separated by a time interval of ~2 s. 
The resulting sequence, which is used for further G-band radiative 
diagnostics, covers approximately 30 minutes of physical time, cor- 
responding to ~4 — 6 granular lifetimes. G-band images were pro- 
du ced using t he radi ative diagnostics technique briefly described 
bv lJess et alj ( 1201 2al) . The plasma velocity field corresponding to 
these images was measured at a height of about 500 km above the 
continuum formation layer (see Figure[2}. 



3 ANALYSIS AND RESULTS 

iKevsetalJhouh employed an automated detection and tracking 
algorithm to estimate the velocities of 6000 G-band MBPs. How- 
ever, an automated detection and tracking algorithm is extremely 
hard to employ for the Ca II K data set. Due to the reverse granu- 
lation inherent in the Ca II K images, the application of intensity 
thresholding for the identification of MBPs proved unsatisfactory. 
The velocities of Ca II K MBPs were instead de termined using 
the L ocal Correlation Tracking (LCT) technique dMatloch et al. I 
I2OIOI) . This was utilised along with manual tracking to ensure 
accurate velocity estimates. As shown in the G-band analysis of 
iKevs etalj|20T ij), the LCT technique produces very similar results 
to the automated tracking algorithm. 

The ambiguity associated with their boundaries did not allow 
a one-to-one correspondence between all the MBPs in the G-band 
and in Ca II K. Approximately 100 MBPs were analysed in Ca II K 
with half of them having unambiguous G-band counterparts. The 
velocity and area of each MBP was subsequently evaluated. As 
mentioned previously, the actual value for the area of an MBP in 
Ca II K can be difficult to determine due to the ambiguity of the 
MBP boundaries. The error in the determination of MBP surface 
area was established through the calculation of mean absolute mea- 
surement errors, and found to be approximately ±23%. Difficulties 
in defining an accurate area for Ca II K MBPs can be clearly seen 
in Figure[T](b). 

Our sample of 100 Ca II K MBPs revealed an average velocity 
of 1.4 km s^^, with an estimated average area of 52 000 km^ The 
area of Ca II K MBPs follows a log-normal distribution similar to 
the G-band (Crockett et al. 2010), with the Ca II K area distribu- 
tion shifted to higher values with a broader peak (Figure [3]a). Ve- 
locity distributions for the observed and simulated G-band datasets 
and Ca II K are shown in Figure [3] (b). The transverse velocities 
of Ca II K MBPs are enhanced compared to the G-band (an aver- 
age of 1 kms^^ for 6236 MBPs), with the mean absolute measure- 
ment error in both velocity values estimated to be 0.5 km s^^. Their 
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Figure 2. Sample images from the simulated data sets, (a) Snapshot of a 12x 12 Mm G-band simulation, (b) The con'esponding velocity domain 500 km 
above the continuum formation layer for the image shown in (a). Velocity scale is in kms^^. A sample, isolated MBP, is contoured red in (a) and (b). (c) 
shows a set of expanded panels depicting the evolution of a MBP found in the white box in (a). The top row shows the G-band images of the MBP advancing 
in time by ~44 s between each panel from left to right. The bottom row shows the corresponding velocity domain images for the MBP. The final panel shows 
the MBP as seen in (a) and (b). 



size also appear to be larger by about 63%. The distributions were 
analysed in more detail by selecting 50 Ca II K MBPs and their 
corresponding G-band counterparts. Again, the average transverse 
velocities are 1.8 kms~^ and 0.8 kms~^ in the Ca II K and G-band, 
respectively. The corresponding average sizes are 51 800 km^ and 
31 300 km^. This shows a 65% increase in the size of the MBPs 
as they extend upwards, in agreement with ljess et alj ( l2010a) . who 
found an expansion of ~70% in their study of MBPs manifesting in 
the Na I Di line core. Of these 50 MBPs, ~80% are isolated and the 
remiander are groups of MBPs. On closer inspection the isolated 
MBPs have higher average velocities in both G-band and Ca II K 
(0.8 kms^^ and 1.9 kms~^ respectively) compared to the groups 
of MBPs (0.6 km and 0.3 km s"^ in G-band and Ca II K respec- 
tively). As is expected, MBPs found in a group have a larger aver- 
age area in both G-band and Ca II K compared to isolated MBPs. 
Also, the change in area between G-band and Ca II K is more sig- 
nificant for isolated MBPs (~57% change for MBP groups whereas 
nearly twice as large on average for isolated MBPs.) This is likely 
th e result of a higher field strength found in groups of MBPs as seen 
in ljess et al. (2010a). The corresponding velocities for a sample of 
10 of these MBPs are given in Table[T] 



We have also studied the MBPs visible in the Na I Di line. 
The scanning capability of IBIS allows us to study the MBPs both 
in the core and wings of the line, albeit at a reduced temporal res- 
olution of 39.7 s. This resulted in a smaller sample of 50 MBPs 
present in both the wing and line core (Figure [T](c) and (d)). The 
average velocity for the core of the Na I Di line was found to be 
only 0.7 km s^^ and the average for the 486 MBPs analysed at the 
wing of the line 0.7 kms~^, with a maximum value of 6 kms~^. 
This maximum value is close to that found for the G-band. How- 
ever, there is a slight discrepancy between the other velocities es- 
tablished for Ca II K, G-band and Na I Di. This is probably a result 
of the combination of a smaller sample size and the lower temporal 
resolution of the Na I Di data sets. The errors in the velocity values 
in the core and the wing of the Na I Di line are estimated to be 
0.6 km s^^ and 0.7 km s^^ respectively. 

The observational results were compared with radiative MHD 
simulations pro duced with the MU RaM code. These simulations 
are described in |jessetal.l ( [2012bl) and have a cadence of ~2.l s. 
We have 724 G-band MBPs in these 200 G simulated images, 
which have an average velocity of 1.3 km s^^, a maximum velocity 
of 6.9 kms^^, and a mean lifetime of 110 s. The simulated MBP 
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Table 1. Summary of MBP characteristics directly compared at various lieights. 



MBP 


Ca II K 


Ca II K 


G-band 


G-band 


No. 


Vel.(kms 1) 


Area(km^) 


Vel.(kms-i) 


Area(kin^ ) 


1 


0.74 


95,000 


0.75 


35 000 


2 


3.26 


66,250 


1.10 


40 000 


3 


2.09 


88,750 


1.62 


40 000 


4 


1.27 


131,250 


0.46 


35 000 


5 


1.36 


88,750 


0.81 


30 000 


6 


1.19 


87,500 


0.23 


50 000 


7 


2.04 


142,500 


0.13 


102 500 


8 


0.82 


113,750 


0.51 


52 500 


9 


2.77 


116,250 


2.54 


112 500 


10 


0.26 


156,250 


0.68 


110000 



Table 2. Summary of MBP characteristics for observations (ROSA/IBIS) and simulations (MURaM). 



Data Set 


FOV 
Size 


Av. 
Cadence(s) 


No. 
MBPs 


Av. Vel. 

(kms-i) 


Max. Vel. 

(kms-i) 


Av. Area 
(km2) 


G-band 

(ROSA) 


70" X 70" 


2.1 


6236 


1 


7 


32 000 


Ca II K 

(ROSA) 


70" X 70" 


4.2 


100 


1.4 


10.2 


52 000 


NalDi 
(Wing-/B/5) 


83" X 83" 


39.7 


486 


0.7 


6 


51200 


NalDi 
(Core-ffi/5) 


83" X 83" 


39.7 


50 


0.7 


2.4 


73 500 


G-band 
(MURaM) 


16" X 16" 


2.1 


724 


1.3 


6.9 


26100 


Velocity Domain 


16" X 16" 


2.1 


724 


1.8 


15 


N/A 



(MURaM) 



velocities have an error of approximately 0.5 kms .A summary 
of results from the observed and simulated data sets is given in Ta- 
blelH 

The MURaM simulations also allow us to study the corre- 
sponding velocity field at a height of 500 km above the contin- 
uum formation layer. Transverse velocities could therefore be es- 
tablished for each point in the data set along with an inclination 
angle. The velocity maps were compared to the MBPs detected in 
the simulated G-band images and used to determine an average ve- 
locity for MBPs at this height of 1.8 km s^^, with a maximum value 
of 15 kms^^. 

Our simulations reproduce the measured G-band and Ca 11 K 
values, albeit with slightly elevated velocities when compared to 
the observations. However, if we were to consider only the 50 un- 
ambiguous MBPs identified in both bandpasses, we can see that 
the velocities for MBPs found higher in the atmosphere are the 
same. The simulated images are from a format i on hei ght of ap- 
proximately 500 km whereas iBeebe & JohnsonI ( 1 19691) state that 
Ca II K is formed at a height less than 1300 km. Also, the pres- 
ence of reverse granulation (as displayed in Figure [TJ would sug- 
gest that the dominant contribution to the Ca II K filter com es from 
the up per photosphere/lower chromosphere. A review by iRuttenl 
( l2007h highlights the complexity of determining line formation re- 



gions from narrow band pass filters and suggests that Ca II K is 
formed at the lower chromospheric level. Although the simulated 
images do not have a formation height which is explicitly the same 
as the observed Ca II K core images, we think that a comparison 
can still be made. The relatively broad bandpass of the Ca II K filter 
make it difficult to assign a specific height of formation. However, 
there is a slight discrepancy between the values found at the forma- 
tion height of the G-band, with the velocities of the observed data 
having an average of 0.7 kms~^ compared to 1.3 kms~^ in the 
simulated images. This is possibly due to the far superior sample 
size used to ascertain this value in the simulated images. Both the 
simulated and observed datasets show that the average velocities of 
MBPs increase with height. This is probably a direct consequence 
of the conservation of momentum co mbined with a de crease in the 
density. However, it has been shown jjess et al]l2010al) that the ex- 
pansion of flux tubes as they extend upwards is affected by the field 
strength of the tube. Higher field strengths lead to less expansion of 
the flux tube at higher levels. This could also affect the velocities 
of the MBPs, but more studies would be required to ascertain what 
effect the field strength has, if any, on the estimated velocities. 

To understand the nature of these flux tubes, and how they 
transfer energy to higher levels, we look for changes in the ve- 
locities of MBPs between different levels as a function of time, 
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Area (km^) Velocity (km/s) 

Figure 3. The area and velocity distributions of the observed and simulated datasets. (a) The normalised area distributions of the data sets; the black line is the 
observed G-band, the red line is the simulated G-band and the blue line is the Ca II K data. The log-normal distribution of the two G-band data sets is clear. 
However, the smaller number of data points in Ca II K means that this distribution is not as clear. The dashed-blue line shows a possible ideal log-normal 
curve distribution for Ca II K given better count statistics. It is evident that the Ca II K MBPs have larger areas than their G-band counterparts, (b) displays a 
histogram of the velocity distributions for observed G-band {black), simulated G-band {red) and Ca II K {blue). The G-band distributions are fairly similar and 
that for Ca II K illustrates elevated MBP velocities. 



which may indicate how rigid these flux tubes are. This study was 
first performed on the simulations as they contained the velocity 
information for higher layers, making the analysis easier. From this 
analysis an average delay of ~4 s was found between the G-band 
simulations and a height of 500 km above the continuum formation 
layer where the velocity domain is taken. This may indicate that the 
flux tubes between the layers are not rigid. However, as the spread 
of values for the delays in velocity ranges from ~2 — 8 s it suggests 
that the rigidity of these magnetic flux tubes can vary. Both the field 
strength and geometry of the individual tube can play a role in its 
flexibility. 

The average magnetic field strength of G-band MBPs in the 
simulations at the continuum formation level is about 950 G, while 
that of the corresponding MBPs detected at 500 km above the con- 
tinuum formation layer is 450 G. For the MBPs studied, the mag- 
netic field strength is therefore reduced by ~50% between the G- 
band images and the height in the atmosphere where we take the 
velocity measurements. Over the same range the density is reduced 
by approximately a factor of ten. 

The short time delays in the simulated images at the contin- 
uum formation layer and a height of 500 km would imply that en- 
ergy between the two layers is transported at a very high rate. Given 
that the sound and Alfven speeds are in the range 11 — 12 km s^^ 
and 10 — 15 kms~^, respectively, between the two layers, it seems 
unlikely that the observed delay is due to the movement of the foot- 
points in the simulation. A more likely scenario may be that the 
delays are produced by oblique granular shock waves, as reported 
bylvkas et al. (2011). 

A similar analysis was performed on the observed G-band and 
Ca II K datasets. The observations show an average delay of 29 s in 
the response of the velocity changes between the bandpasses, with 
the actual values ranging between 21s and 42 s. 



Due to the broad bandpass of the Ca II K filter, it encompasses 
radiation from both the chromospheric emission core and photo- 
spheric absorption profile. It is therefore extremely difficult to as- 
sign a specific height from which this radiation originates. How- 
ever, the reverse granulation present in the Ca II K images indicates 
that emission from the upper photosphere may dominate this band- 
pass. A rough estimate for the distance between the G-band and 
Ca II K emission region is around 500 km. Given a delay of 29 s, 
this suggests a velocity of propagation of ~17 km s^^ between the 
two positions for the flux tubes. This would suggest that the prop- 
agation velocity is close to the sound speed, given typical values 
for the two positions. As the propagation speeds are slower in the 
observations, it suggests that the physical reason behind the delay 
may not be the result of granular shock waves. It is possible that 
the delay between the G-band and Ca II K observations is the result 
of two interwoven factors; the decrease in density and the expan- 
sion of the flux tube with height. The first factor is easily explained. 
With a decrease in the density of the surrounding material, the flux 
tube is less confined and therefore less responsive to the granular 
movements at the photospheric level, thus resulting in delays of ve- 
locity shifts. The second factor, which as stated previously is most 
likely a result of the decrease in density with height, probably has 
an added effect on the flexibility of the tubes. As the tube expands in 
the atmosphere, the magnetic flux contained within is less densely 
held in the tube and as a result it is less rigid. This may explain the 
difference in delays seen in the ROSA images. Unfortunately, we do 
not have magnetic field information for both G-band and Ca II K to 
assess how the field strength varies with height in the observations, 
and to investigate if there is any correlation between this and the 
delays in velocity shifts. 
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4 CONCLUSIONS 

We study the properties of MBPs as a function of lieiglit in tlie 
lower solar atmosphere using high resolution multi-wavelength 
imaging and radiative MHD simulations. The results show that 
the transverse velocity and area coverage of MBPs increases with 
height. Velocities of 1 kms~^ are observed for MBPs found in the 
G-band images, with the corresponding Ca II K images yeilding an 
average velocity of 1 .4 km s ~ ^ . The areas of the MBPs show an av- 
erage increase of approximately 63% from the photosphere to the 
chromosphere. 

Radiative MHD simulations reveal similar findings but the 
theoretical values are slightly elevated when compared to the ob- 
servations. Due to difficulties associated with simulating the radia- 
tive emission in Ca II K, we are unable to calculate changes in the 
area and velocity of MBPs in simulated Ca II K images. Instead we 
use the simulated horizontal velocity field and find that the delay 
in the changes of velocity bewteen the two layers is approximately 
4 s. We speculate that this delay may be attributed to horizontally- 
propagating oblique granular shock waves. 

The observations indicate an average delay of 29 s between 
the G-band and Ca II K. We surmise that the difference between the 
observations and the simulations is due to the height of formation 
of Ca II K, although the relatively broad bandpass of the Ca II K 
filter makes this height difficult to ascertain. It should also be noted 
that the difference in properties of MBPs in the observations and 
the simulations could be due to the different nature of turbulence 
conditions in the Sun and in the simulations, where the Reynolds 
number is very low. This is likely to affect the interaction of the 
flux tubes and convection at different altitudes. We have presented 
one of the larger-scale studies of the properties of MBPs at multiple 
heights. The results show that the transverse velocity of MBPs in- 
crease with height, most likely as a result of the decreasing density 
at higher levels and the conservation of momentum. 
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